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Phases: geometric; dynamic or topological 

Mechanical effects of light on atoms, molecules, and ions 

Matter waves 



We propose a semi-classical interpretation of the geometric scalar and vector poten- 
tials that arise due to Berry's phase when an atom moves slowly in a light field. Starting from the 
full quantum Hamiltonian, we turn to a classical description of the atomic centre-of-mass motion 
while still treating the internal degrees of freedom as quantum variables. We show that the scalar 
potential can be identified as the kinetic energy of an atomic micro-motion caused by quantum 
fiuctuations of the radiative force, and that the Lorentz-type force appears as a result of the 
motion-induced perturbation of the internal atomic state. For a specific configuration involving 
two counter-propagating Gaussian laser beams, we relate the geometric forces to the radiation 
pressure and dipole forces known from quantum optics. The simple physical pictures provided by 
the present analysis may help for the design and the implementation of novel geometric forces. 
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Cold atomic gases are considered as efficient simulators 
of quantum condensed matter systems (for a review, see 
e.g. [1]). The confinement potential and the atomic in- 
teractions can be tailored almost at will, allowing for ex- 
ample to mimic with atomic vapours situations encoun- 
tered for electrons in solid state materials. A major step 
ifi the implementation of these simulators is the possibil- 
ity to apply a gauge field to the cold atomic gas in or- 
(Jer to model the vector potential appearing when charged 
particles are placed in a magnetic field. Up to now such 
gauge fields have been mostly obtained by rotating the 
gas [2]. In this case the transformation to the rotating 
frame indeed corresponds to giving the particles a ficti- 
tious charge, and applying an effective uniform magnetic 
field. Another method consists in using so-called geomet- 
ric potentials which can considerably extend the range of 
gauge fields realisable in neutral gases. In particular, they 
offer the possibility to produce non homogeneous or time 
dependent effective orbital magnetism. 

Geometric potentials [3-6] generally result from Berry's 
phase [7] that appears when particles with an internal 
structure move slowly enough that their fast internal dy- 
namics adiabatically adjusts to the centre-of-mass motion. 
In quantum optics such geometric potentials can be gener- 
ated using laser beams to split the atomic internal energy 



levels. This was suggested in [8,9], and a first experimen- 
tal investigation was presented in [10]. The concept of 
geometric potentials can be extended to simulate gauge 
fields that arc more elaborate than the C/(l) symmetry of 
electromagnctism. Using appropriate laser configurations, 
one can in principle implement the general ideas outlined 
in [11] to generate non-Abelian gauge fields [12,13]. 

In spite of numerous investigations of possible geomet- 
ric potentials in quantum optics, a simple physical inter- 
pretation of the forces appearing due to the gauge fields 
seems still to be lacking. Of course these forces ultimately 
arise from the exchange of momentum between light and 
atoms, but this generic process can lead to different phys- 
ical mechanisms. Since identifying these mechanisms may 
help to design future configurations, we propose in this let- 
ter a semi-classical analysis of geometric potentials. Our 
approach is directly inspired by the formalism used to cal- 
culate the standard radiative forces acting on an atom 
placed in a laser beam. We first present the general semi- 
classical derivation of the two (scalar and vector) geomet- 
ric potentials acting on an atom. We then discuss the 
various physical pictures that emerge for the concrete im- 
plementation that has been proposed in [14]. 

We start with a brief reminder of the standard formal- 
ism of geometric potentials, in which both internal and 
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external (centre-of-mass) atomic degrees of freedom are 
treated using quantum mechanics. The relevant internal 
dynamics is described in an A^-dimcnsional Hilbcrt space. 
The atom can be submitted to static electric or magnetic 
fields as well as to time-dependent electromagnetic fields. 
The fields are supposed to be in a coherent state so that 
they can be described by classical functions. Assuming 
that the time-dcpendcncc of the atom-field interaction can 
be eliminated using the rotating-wavc approximation, the 
eigenstates of the atom-field coupling form an orthogonal 
basis {\ipj{r)), j = 1, . . . , N} of the internal Hilbcrt space 
of the atom at any point r. We denote Ej{r) the cor- 
responding energies. The Hamiltonian of the problem is 
thus 



TV 

V{r) = Y.^,{r)Q,{r). 



(1) 
(2) 



Here, M is the atomic mass, p ^ —iTi'V the centre-of-mass 
momentum operator and Qj{r) = \ipj{r)){ipj{r)\ the pro- 
jector onto the j-th internal eigenstate. Suppose now that 
the energy of one of the internal eigenstates, say Itpi), is 
well separated from the other ones. We choose the initial 
internal atomic state equal to \'ipi{r)) at any point r and 
suppose that the atom moves slowly enough for the adi- 
abatic theorem to hold. The internal state then remains 
equal to {ipi) at any time, and the energy Ei{r) plays 
the role of a potential energy for the centre-of-mass mo- 
tion. In addition, the geometric phase accumulated by the 
atom gives rise to additional vector A{r) and scalar U{r) 
potentials so that the atom Hamiltonian in the adiabatic 
approximation [3-6] reads 



H.. 



adiab. 



ip-Mr)y 

2M 



+ Ei{r) + U{r) , 



with 



A{r) 
U{r) 



i?i(V'i(r)|V^i(r)) , 



2M 



Y^\{Mr)\V^,{r))t 



(3) 

(4) 
(5) 



i/i 



where we note by convention |V'0(r)) = V {\il:{r))). 

The goal of this letter is to provide a simple physical 
interpretation of these potentials within the framework of 
a semi- classical analysis. Here, the term semi-classical 
refers to the fact that we describe the atomic internal de- 
grees of freedom using quantum mechanics, but we treat 
classically the centre-of-mass motion. Within this approx- 
imation we want to recover the equation of motion 



fill 

M^ = -VE,{r) 
at 



VUir) +v X B{r) 



(6) 



corresponding to the Hamiltonian (3). This equation of 
motion involves three forces. The first one is simply the 



gradient of the energy Ei of the occupied level. The sec- 
ond one originates from the scalar potential [/(r), and the 
third one has the structure of a Lorentz force for a charge 
g = 1 in an effective magnetic field B{r) = V x A{r). 

The main tool for the semi-classical analysis is the force 
operator acting on the iV-dimensional Hilbert space: 



F{r) = -VV ^-J2 (^(Ej) Qj + Ej V{Qj) 



(7) 



Knowing the internal state |0) for an atom at point r (a 
concept which is valid within the semi-classical approach) , 
wc will be able to calculate the average force {(l)\F{r)\(f>) 
and the correlation functions of the force operator. This 
semi-classical approach has been very fruitful for the study 
of the radiative forces acting on an atom irradiated by 
laser beams [15], and its connection with a full quantum 
description of the atomic motion in laser light is well es- 
tablished [16]. 

Origin of the scalar potential. In this section 

we consider an atom with zero centre-of-mass velocity, so 
that the Lorentz force in (6) is also zero. The atom internal 
state is supposed to be [(/)) = [i/ji) and we immediately get 
from (7) 

(F{r)) = -VE, , (8) 

where we used that the states {\i^j{r))} form a normalised 
orthogonal basis at any point r. We thus recover the first 
term in eq. (6) but not the force originating from the scalar 
potential U{r). This could be expected since the expres- 
sion (5) of the scalar potential involves the atomic mass 
M, which does not enter in the expression (7) of the force 
operator F. The scalar potential can still be recovered 
within the semi-classical approach, as we show now, pro- 
vided one goes one step beyond the mere calculation of 
the average force. 

The starting point of our reasoning consists in noting 
that the state [^/ii) occupied by the atom is an eigenstate 
of the coupling Hamiltonian V, but not an eigenstate of 
the force operator F. Hence, (.F^) 7^ {F)'^ or, in physi- 
cal terms, the force acting on the atom undergoes quan- 
tum fluctuations around its average value. As we will see, 
these fluctuations occur at the typical Bohr frequencies 
LOji = {Ej — Ei)/h of the internal atomic motion. For the 
adiabatic approximation to hold, these frequencies have 
to be much larger than the characteristic frequencies of 
the external atomic motion. Consequently, the quantum 
fluctuations of the force cause a fast micro-motion of the 
atom, similar to the one arising for charged particles in 
a Paul trap [17]. The kinetic energy of the micro-motion 
then plays the role of a potential for the slow motion of 
the atomic centre-of-mass [18]. We demonstrate below 
that this kinetic energy coincides with the scalar potential 
C/(r). 

The fluctuations of the force operator F are charac- 
terised by the symmetrised correlation function of the op- 
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erator 6F ^ F — (F) , calculated in the Heisenberg picture: 
C{t, t') = i {5F{t) ■ SF{t') + 5F{t') ■ dF{t)) . (9) 

Since the average is taken in an eigenstate of the Hamilto- 
nian, C{t, t') depends only on the time difference t = t — t' 
and we obtain after some algebra 






with 



(10) 



(11) 



In order to understand the consequences of these fluctua- 
tions, consider a classical particle submitted to a stochas- 
tic force F{t) such that F(t) = and F(t) ■ F{t') = 
C{t — t'). The Fourier transform f{oj) of F{t) is such 
that f{uj) = and f*{uj) ■ f{uj') = 5{uj - u') B{uj), where 



B{u) 



'U 



Cj {Siuj + Uji) + S{uj - Uji)) (12) 



is the Fourier transform of C(t) . The solution of the equa- 
tion of motion p = f is 



pm^l 



/(") 



duj 



(13) 



and has a zero average. However, the average kinetic en- 
ergy is strictly positive and equal to 



2M 



2Afw2 



duj . 



(14) 



In the explicit calculation of (14), the contribution of the 
w~2 denominator cancels out the transition frequencies 
LOji that appear in the expression (11) of Cj. Finally 
one exactly recovers the result (5) for the scalar poten- 
tial. This validates the interpretation of this potential in 
terms of the kinetic energy of the atomic micro-motion. 

The above result sheds new light on the Hamiltonian 
(3) of the full quantum description. We can now interpret 
the term Mv'^/2, with Mi) — p— A, as the kinetic energy 
of the slow centre-of-mass motion, whereas the kinetic en- 
ergy of the fast micro-motion builds the scalar potential 
U{r). It is also interesting to connect the present anal- 
ysis of the scalar geometric potential with the intriguing 
problem of a two-level atom moving around the node of 
a standing wave. In the latter case, it is found that the 
average force acting on the atom is zero, as expected since 
the light intensity vanishes at the nodes. However, the 
atomic momentum diffusion coefficient, which is also re- 
lated to the correlation function of the force operator, is 
non zero [15, 19]. 

Origin of the Lorentz force. We now consider 

the case of a slowly moving atom and calculate the av- 
erage of the force operator (7) at first order in velocity. 



More precisely, we assume an atom initially at rest in the 
internal state 10) = |^i), that is adiabatically set in mo- 
tion to reach a velocity v. Because of this motion the in- 
ternal atomic state contains some admixture of the other 
cigenstatcs lijjj) and the average force is different from the 
zero- velocity result. 

We write the internal state as |0) —J2i '^j iV'j) a-nd solve 
the Schrodinger equation as a power scries in velocity. The 
procedure detailed in the appendix gives the coefficients 
a, at first order in v: 



aj (i) ~ ih 



MV^jIWi) ,B,i 



/h 



E, - E, 



U + 1) • (15) 



The calculation of the average force at first order in v is 
also outlined in the appendix and leads to 



(F) = i?i(VV'i| {v ■ IVV-i)) + c.c. 



(16) 



One can readily check that this expression coincides with 
the Lorentz force v x B appearing in eq. (6). 

This way of recovering the Lorentz force is very rem- 
iniscent of the general semi-classical calculation of the 
velocity-dependent radiative forces [15,16]. There is one 
important difference, however. In the latter case one gen- 
erally finds F -v 7^ 0, which is at the origin of laser cooling 
(for example via the Doppler or Sisyphus mechanisms). 
In the particular case considered here, no photon sponta- 
neous emission process occurs, dissipation is absent, and 
we are left with a Lorentz force of geometric origin. 

Illustration for a particular atom-laser configura- 
tion. — We now relate the geometric forces in quantum 
optics to the known radiative forces that generally act on 
an atom irradiated by one or several laser beams. We 
thus explain the mechanisms at the origin of the geomet- 
ric potentials in terms of exchange of momentum between 
atoms and light. For this purpose, we turn to the specific 
configuration sketched in fig. 1, which has been proposed 
in [14]. A three-level atom with two degenerate ground 
states \g±) and an excited state |e) is irradiated by two 
counter-propagating laser beams. The beam propagating 
in the +y (resp. —y) direction drives the transition g^ <-> e 
(resp. g- <-^ e). The two beams have the same waist if 
and the same intensity. They arc offset with respect to 
the y axis by a distance ±6/2, where b is typically of the 
order of w. 

Using the rotating-wave approximation the atom-light 
coupling can be written as 

Vir) = -?iA|e)(e| + ^ {h^,{r)\e){g,\ +h.c.) . (17) 



The Rabi frequencies k± are given by 



(18) 



where k is the modulus of the wave vector of the laser 
beams. We assume that the Raman resonance condition is 
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Fig. 1: Atomic level scheme and laser configuration proposed 
in [14] to generate a vector and a scalar potential. 



satisfied between the two ground levels g±, and we denote 
A the detuning of the laser frequency with respeet to the 
(7-1- <-> e transition. 

A well known eharacteristies of this configuration is that 
one of the three eigenstates of V{r) involves only the two 
ground states \g±) and has zero overlap with the excited 
state |e) [20-22]. An atom prepared in this so-called non- 
coupled or dark state will not undergo any spontaneous 
emission process, which is a key feature for the practical 
use of geometric potentials. The non-coupled state is 



\^i 



n 



where we set fl 



n 

|2U/2 



■|5-) , 



(19) 



-I- ]k-|-] ) ' . The corresponding 



energy at any point is _Ei(r) = so that the first term on 
the right hand side of (6) vanishes. The scalar potential 
and the effective magnetic field are [14]: 



ui.) ^ '^^^iP^ci.), 



B{x) = 



2M 
5- G{x) Uz , 



(20) 
(21) 



where G(x) = coshr"^ {2xb / ■uP') . 

We start our analysis with the scalar potential U{r). 
When X -^ ±00, it tends to zero as expected from its inter- 
pretation as the energy of a micro-motion, which vanishes 
when the atom sits outside the laser beams. The value of 
U at r = can also be recovered easily. At this point the 
two components of the force operator are 



F. = -Fo (le)(.ga| + |.9a)(e|) , (22) 

(23) 



Fy = -^Fi(le)(5a|-|5a)(e|) 



where Fq ~ \/2fiKb/w'^, Fi = y2fikK and where we 
have introduced the antisymmetric combination \ga) = 
i\g+) - \g-))/V2. At r = 0, J^i) = \ga) so that (F^^y) = 0, 
which directly follows from Ei{r) ~ 0. The correlation 
function of the force is readily calculated; choosing A = 
for simplicity, we find C{t) = {Fg + F^) cos{V2 kt). The 
kinetic energy of the micro-motion induced by this oscil- 
lating force is {Fq + Fl)/{AMiP) and indeed coincides 
with the general result (20) at x = 0. For realistic laser 
configurations the waist w and the offset h are large com- 
pared to fc~^ so that the dominating contribution is due 



to Fi , corresponding to a micro- motion directed along the 
propagation axis y of the beams. 

The two operators F^ and Fy given in (22-23) are the 
so-called dipole force operator and radiation pressure force 
operator, respectively. The dipole force F^ originates from 
the intensity gradient of the laser beams along the x axis. 
In terms of photon momentum exchange, it can be under- 
stood as a redistribution of photons between the various 
plane waves that contribute to the formation of the inten- 
sity gradient (see e.g. [19]). The radiation pressure force 
Fy originates from the phase variation e^**^^ of the laser 
beams along the y-axis. It leads to changes of atomic mo- 
mentum by ±hk when the atom absorbs a photon from 
one of the laser beams. Interestingly, the cancellation of 
the average force acting on an atom at rest at r = in 
the internal state \ga) results from a destructive interfer- 
ence. For example, (Fx) = because the state \ga) is a 
combination with equal weights of the two eigenvectors 
\x±) = (Iffo) ± |e))/>/2 of the force operator F^. The 
dipole force felt by an atom in state \x±) is TFq, so that 
the average force vanishes for an atom prepared in \ga)- 

We now turn to the case of a moving atom and to the 
interpretation of the Lorentz force. As the force is linear 
in velocity, we can study separately the cases of a motion 
parallel to the x and to the y axis. We first consider an 
atom moving along the y axis with velocity Vy . Its internal 
state differs from the value \ga) of an atom at rest, and its 
expression at first order in Vy, as deduced from (15), is 



1.9a; 



y I 
V2k 



(24) 



where we have again chosen A = for simplicity. This 
state can also be calculated directly by switching to the 
atom rest frame, which amounts to adding the small per- 
turbation SV = kvy{Q- —Q+) to the atom-light coupling, 
where Q± is the projector on the state |<?±). One can 
check that the state \(j)) is an eigenstate of the perturbed 
coupling V + dV with the same eigenvalue 0. The average 
force acting on an atom in the state (24) is not zero any- 
more. The perfect balance between the two dipole forces 
±Fo is now broken to the benefit of —Fq (-l-i^o) if Wj, > 
{vy < 0). The explicit calculation of {cf>\Fx\(j>) is immedi- 
ate and yields exactly the Lorentz force value —2'hkbvy/w'^ 
obtained from eq. (21). In the case of a motion along the 
y axis, the Lorentz force is therefore a direct consequence 
of the dipole potential. Remarkably, the Lorentz force is 
independent of the Rabi frequency k although the dipole 
force amplitude Fq is proportional to k. The reason is 
that the "rotation" angle of the state \(j)) with respect to 
the non-coupled state scales as k~^ [see eq. (24)] and this 
scaling exactly compensates for the k dependence of Fq. 

We finally consider an atomic motion along the x axis, 
with a Lorentz force along the y-axis, originating thus from 
the radiation pressure force operator Fy. This case could 
be analysed along the lines of eqs. (15-16), but it is actu- 
ally more instructive to look at the cumulative effect of the 
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Lorentz force, i.e. the change of momentum along y, when 
the atom crosses the two laser beams. Suppose that the 
atom is initially located at time ti at a point xi < with 
\xi\ ^ w. At this location k+/k^ — exp{2xib/w^) ^ 1, 
so that eq. (19) leads to \ipi) ~ \g+). At a later time t2 the 
atom has reached the point X2 ~ —xi, where Jt/ii) ~ Iff-). 
The momentum change caused by the Lorentz force is 
aligned with the y axis and its value is 

rt2 /•X2 

Apy = — / BzVxdt = — I Bz{x)dx 

Jti Jxi 

= hk[tSinh{2x2b/w'^)-i&nh{2xib/w'^)\ .(25) 

With the above assumption that |a;i| = \x2\ 2> w and 
b ^ w, we obtain 

Apy ~ 2hk . (26) 

This expression has a clear physical interpretation in terms 
of photon absorption and stimulated emission. When the 
atom moves from xi to X2, its internal state rotates adia- 
batically from Igj^.) to \g-)- This happens by the absorp- 
tion of a photon from the wave propagating in the +y di- 
rection, driving the ff-i- — > e transition, and the stimulated 
emission of a photon in the wave propagating in the —y 
direction, driving the e ^ g- transition. The origin of the 
Lorentz force in this case is therefore closely related to the 
so-called STIRAP process (Stimulated Raman Adiabatic 
Passage) [23]. 

Discussion and conclusion. We have given here 
a semi-classical interpretation of the geometric forces that 
act on a slowly moving particle with multiple, spatially 
varying, internal energy levels. Such a situation is encoun- 
tered in quantum optics when an atom moves in the light 
field created by several quasi-resonant laser beams. The 
main assumption is that the atom follows adiabatically one 
of its internal energy levels which is supposed to be well 
separated from all the other ones. Under these conditions 
the atomic centre-of-mass motion involves a scalar and a 
vector potential. The scalar potential originates from the 
kinetic energy of the micro-motion of the atom, under the 
action of the quantum fluctuations of the radiative force 
operator. The Lorentz-type force associated with the vec- 
tor potential results from the perturbation of the atomic 
internal state due to the slow atomic motion. This changes 
the expectation value of the force with respect to an atom 
at rest, but does not induce any dissipation. The way we 
obtain our results is very reminiscent of the general deriva- 
tion of the radiative forces created by quasi-resonant laser 
beams. Taking as an example the configuration proposed 
in [14], we have related the intervening geometric forces 
with either the radiation pressure force or the dipole force. 

Our approach can be viewed as a quantum version (for 
the internal degrees of freedom) of the fully classical re- 
sults of [24], where the authors analysed the motion of a 
particle with a permanent magnetic moment placed in a 
strong, non homogeneous magnetic field. There, a rapid 
oscillation of the magnetic moment superimposed onto a 



slow secular motion was found. This rapid oscillation gives 
rise to a scalar potential similar to the one of interest 
here. When the particle with its magnetic moment ini- 
tially aligned with the magnetic field was set in motion, it 
was found that the magnetic moment acquires a non-zero 
component in a direction perpendicular to the magnetic 
field. The Lorentz-type force that emerged in this situa- 
tion is also very similar to the present one. 

To conclude, geometric potentials may constitute in the 
future a key ingredient for the realisation of a general 
'quantum gas toolbox' that allows one to address various 
open problems of many-body physics. We hope that the 
simple physical pictures provided by the present analysis 
may help for the design and the implementation of novel 
geometric forces in this context. 

Appendix: Derivation of the Lorentz force. In 

this appendix we detail the calculation of the perturbed 
state [(/)) of a slowly moving atom and the resulting force 
{F). At time the atom is at rest in the internal state 
10) = IV'i)- We suppose that it is uniformly accelerated 
during a time T to a velocity v. v{t) ^ v t/T for < i < 
T . Our aim is to calculate \(j)) and {F) at time T. 

The general expression of the internal atomic state is 
10(f)) = ^jOij{t)\'ijjj{r{t))), leading to the average of the 
force operator (7): 

(F) = -^Vi?, |a,f 
j 

+ J2iEk - E,) ala, {^k\V4',) • (27) 

To calculate the force at first order in v, we need all coef- 
ficients ttj also at first order. Using the Schrodinger equa- 
tion i?i|(/)) = T/(r(t))|0) and 10) = Ej "jlV'j)+aj ■"■[VV'j), 
we get the corresponding equations of motion: 

dj = -lEjaj/h + ^akv- {ipj\V^k) ■ (28) 

k 

At order zero all a^'s are zero except ai(t) = 
ex-p{—iEit/h). At order one we obtain for j ^ 1 

a,{T) = -V ■ (^jIV^i) e-^^^^/^ ^ ^^{E,-E,)t/n 1 ^^ ^ 

(29) 
where [V'j), \'4'i)i Ej and Ei are taken at order zero in v, 
hence at the location r of the atom at time T. Assuming 
that the atom is adiabatically set into motion, i.e. T{Ej — 
Ei)/h-» 1, we get 



«j(r) 



^^MMZM ^-^E,T/n 



(30) 



E, - E, 

At order one in v the equation of motion for ai is 

di = -i{Ei — V ■ A) ai/fi , (31) 

whose solution is a number of modulus 1 . The two results 
(30) and (31) entail that the first part of {F) in (27) has 
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no first order component in v since the contributions of the 
aj's for j ^ 1 are at least of order 2, and the contribution 
of ai is independent of v. In the second part of (27), the 
only relevant terms are those where one of the two indices 
k or j equals 1. Applying the closure relation and keeping 
the terms linear in velocity, we finally reach the result 
(16) for the average force. We note that our procedure is 
similar to the original derivation of the geometric phase [7] , 
which emerges from the term v ■ A in (31). 
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